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EFFECTIVE SPECIFIC IMPULSE OF EXTERNAL NUCLEAR 
PULSE PROPULSION SYSTEMS 
byThaineW. Reynolds 
Lewis Research Center 

SUMMARY 

A simple self-similar flow model for an external nuclear pulse propulsion scheme 
has been used to investigate factors affecting the effective specific impulse of such a 
system. There are three principal factors that control the effective specific impulse: 

(1) the mean propellant velocity, (2) the fraction of total propellant flow f , which inter- 
cepts the pusher of the vehicle, and (3) a mass loss factor f m , which accounts for other 
mass necessarily ablated from the pusher plate of the vehicle, but which has been as- 
sumed herein to contribute little impulse to the vehicle in doing so. 

Based on the model used, the following conclusions Were drawn: 

1. There is an optimum pulse energy for a given system (i. e. , a given pusher 
diameter, and opacity of ablated pusher material) to yield a maximum specific impulse. 

2. Increasing the mean propellant velocity does not necessarily result in an in- 
creased effective specific impulse for a given system. 

3 

3. Mean opacities of 10 square meters per kilogram or above appear necessary to 
approach the achievement of maximum effective specific impulses. 

4. Increasing the vehicle size (i. e. , increasing pusher diameter) leads to higher 
values of f c but lower values of f m . The resulting effective specific impulse tends to 
increase if the pulse energy is kept at the optimum value. 


INTRODUCTION 

The work on the nuclear pulse propulsion scheme proposed initially about 1955 
(ref. 1) was carried out under the project name Orion (ref. 2). The original concept in- 
volved the use of fission reactions as the energy source. Consequently, the minimum 

size of the energy package associated with this propulsion scheme was of the order of 
12 

4. 18x10 joules (1000 tons - TNT equivalent). The corresponding minimum vehicle 



size to use these pulse energies effectively is quite large. 

There is renewed interest in this mode of propulsion because of the possibility of 
using small fusion microbombs as the energy source. The fusion reaction might be ini- 
tiated in small pellets of fusable material by either an intense laser beam (ref. 3) or an 
intense relativistic electron beam (ref. 4). Nuclear pulse propulsion is of particular in- 
terest because of the stated possibility of achieving high thrust-weight ratios at high ef- 
fective specific impulse values, of the order of 10 000 seconds (ref. 3). 

There is no attempt made in this report to evaluate the feasibility of obtaining the 
fusion reaction by any method. It is assumed that fusion can be attained in any size en- 
ergy release. This report does consider, however, some basic relations for the propel- 
lant flow and heat-transfer problems associated with such a propulsion scheme. A model 
of the process is used to derive some interrelations between the effective specific im- 
pulse and pulse energy size, collimation factor, and mean absorption coefficient of the 
ablated pusher surface material. The general philosophy with which these calculations 
were approached was to attempt to estimate upper limits to the effective specific im- 
pulse. The assumptions made and values of parameters used are discussed in this vein. 


BASIC NUCLEAR PULSE PROPULSION CONCEPT 

A schematic diagram showing the use of externally exploded pulse units to propel a 
space vehicle is shown in sketch (a). A pulse unit containing fusionable material plus 



propellant mass is ejected into position behind the vehicle and the pulse energy triggered. 
The propellant mass expands. A fraction of the propellant intercepts the pusher plate of 
the vehicle and transfers momentum and heat to the vehicle. The heat flux causes some 
ablation of the pusher surface. A succession of such pulses is continued until the de- 
sired total impulse for the mission is obtained. 

In the usual chemical or electric rocket case, all the propellant mass ejected from 
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the vehicle contributes its momentum to the vehicle; thus, the effective specific impulse 
(the impulse per unit weight flow) is a function only of the mean propellant velocity: 



V = v 

Wp g 


( 1 ) 


(All symbols are defined in the appendix. ) 

As is apparent in sketch (a), not all the mass ejected from the vehicle in the nuclear 
pulse case is effective in contributing impulse to the vehicle. For the pulse case, the 
effective specific impulse is 


(I ) = 

sp eff W n + W a 
p a 

where 


( 2 ) 


1(e) = f c MpV (3) 

is the total impulse intercepted by the pusher. The weight of material (other than the 
propellant) which is lost from the vehicle per pulse and which is assumed to contribute 
no effective momentum to the vehicle is represented by W_. In the analysis herein, the 
only such material considered is material ablated from the pusher plate. The mass loss 
factor is defined as 



i M a 
1 + — - 




The effective specific impulse for the pulse system can thus be expressed as 


( 4 ) 


^sp'eff *c ' f m ^sp*base ^ 

where (I s P W = vg. 

To evaluate the specific impulse for this type of system, then, one has to look at the 
mean propellant velocity v that the system can tolerate, the effectiveness with which 
the mass of the pulse can be collimated so as to intercept the pusher plate of the vehicle, 
and the unavoidable mass losses in the system, particularly those lost through interac- 
tion of the high velocity propellant with the pusher surface. 

The analysis that follows considers the characteristics of the expanding propellant 
to derive conditions affecting the interaction between the propellant and pusher surface. 
The principal interactions affecting the performance limitation of the system are the 
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rate of heat transfer leading to pusher ablation, stress limits on the pusher plate mate- 
rial, and pulse unit design as reflected in the maximum amount of collimation that is at- 
tainable. These factors are discussed in the subsequent sections of this report. 


SELF-SIMILAR EXPANSION 

The propellant mass will be at a high temperature and pressure condition immedi- 
ately following the pulse energy release. As this material expands into the vacuum of 
space, some recombinations of ionized and dissociated products will occur. Eventually, 
through expansion, the density drops to the point where interparticle collisions are rela- 
tively unimportant, and the material continues to expand in what is called a "self - 
similar" manner (ref. 5). The characteristic relation of this type of expansion for the 
case herein is 

v(r, t) = - (6) 

t 

In the treatment herein, it is assumed that all the energy absorbed by the propellant 
mass appears finally as kinetic energy With a Maxwellian distribution about the mean ve- 
locity v: 

Vi v 2 (7) 

Also, it is assumed that the condition for collisionless expansion is reached in a 
relatively short distance, compared to the dimensions of the system, so that the expan- 
sion products are effectively emitting from a point source. The adequacy of this as- 
sumption is discussed in references 6 to 8. 

The density at any location and time is given by (ref. 9) 

P< r,t) 5>_ (±f e-fV*) 2 (8) 

tj. 3/2 r 3 \ v/ 

From equation (8), then, the other properties readily follow: 

Mass flow rate per unit area: 

nip(r, t) = p(r, t) * v (9) 


4 



Pressure: 


Energy flow rate per unit area: 


p(r, t) = p(r, tjv 2 


E(r, t) = 

2 

If equation (6) is used, the various property relations obtained are as follows: 


Property 

Equation in terms of 

Equation in terms of E p 

Equa- 

tion 

Density 

M P /t r \ 3 e -(V t ) 2 
, 8 >VW 

2E p /W e -(t r /t) 

7r 3/2 r 3-2\ t ; 

(12) 

Mass 

flow/area 

V /‘r\ 4 e ' (t r /t)2 

//vw 

2E p /‘A 4 e -( t r /t) 2 

WvVV 

(13) 

Pressure 

V 2 M 5 e ' (t r /t)2 

WVJ 

2E P M 5 .-V* 

ww 

(14) 

Energy 

flow/area 

V 3 / t r\ 6 e -V« 2 
2,WW 

V frf/V 4 ) 2 

ww 

(15) 
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Value at maximum 
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E 

in 

= 0. 194 AL. 

m 

3- 


r v 


E 

^m 

= 0.291 -E 
3 
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E v 

E_ 

= 0.241 -E- 

m 

3 
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Some plots of equations (12) to (15) are shown in figure 1 for a particular Ep/jr 
value. The time scales are nondimensional in this figure, the reference time being 
t = r/v. One can note, then, that the time at which the maxima in the various parame- 
ters occur (t ) is a direct function of the separation distance r. Also, the total time of 
pulse interaction is of the order of three times t , so that pulse interaction time also 
varies directly with separation distance. This fact has a bearing on the heat-transfer 
effects as will be seen later. 

Consider the arrival of propellant from a pulse onto a plane normal to the main flow 
direction (sketch (b)). The mass flux at point (z, /3) normal to the plane is 


p(z, 0 ) • v cos 0 = 


_ Mp V COS 0 / v \4 (y/ y)2 


3/2 3 
7 r ' r 
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(b) 


and the pressure on the plane at that point is 


P(z, 0 ) 




cos |3 


_3/2 3 
7 t 1 r 




«2 


,-(v/v) 


(17) 


Typical radial variations of pressure and density are shown in figure 2. The impulse 
per unit area on the plane at (z, 0) in the z- direction is 


K0) 


= 5 > 


V COB fi 


2„ 3 / 2 r 2 


(18) 


The total impulse in the z-direction from flux within a given cone angle 9 is 



But the fraction of the total mass flow which is included within a given cone angle 6 is 

f e =sin 2 ^ (20) 


for isotropic distribution. Thus, considering only the propellant interception factor, the 
effective specific impulse of a system using a circular pusher plate to intercept the mass 
flow within a cone angle 6 is 


(I=J 


sp'eff 


- 1 ( 0 ). 
V 



_v 

g 


( 21 ) 
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If a pulse unit is designed to direct a disproportionate fraction of the total mass into 
a given cone angle, the flow is still assumed uniformly distributed within the cone. The 
parameter C (hereafter referred to as a collimation factor) is defined as the ratio of the 
enhanced total mass in the cone to the amount in the cone if the distribution were isotrop- 
ic. It can be expressed as 



( 22 ) 


With this assumption of the flow distribution, the relations developed for the isotropic 
distribution may be used for collimated flow cases by replacing the propellant mass term 
Mp with the product C Mp. 

If the case with collimated flow is now assumed, the relation for total impulse (eq. 


(19)) is put in terms of pulse energy Ep^= Jr MpV 2 j and C 


1(6) = — cos 

yfn 



s rf\5 


(23) 


The pressure is highest at the center of the plate ()3 = 0) 


P(M) 


2CE p AA5 -(t r /t) 2 

w w 6 


(24) 


At /3 = 0, the maximum pressure at any time is 


CE 

■p m = 0.291 — £ 

o 


(25) 


The effect of two factors (equivalent energy of the pulse and attainable collimation 
factor) involved in the design of pulse unite on the attainable propellant interception fac- 
tor f can be seen by the following development. Combining equations (21) and (23) yields 
c 


^sp^eff _ C 
^spWse 2 yfn 



(26) 


In all subsequent relations where a specific value of maximum pressure is used, a value 

O 

of p m = 6.9x10 newtons per square meter (equivalent to about 100 000 psi) is used. 
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This value represents a reasonable upper limit to the allowable yield stress of materials 
that might be used for the pusher. If this value is used, equation (25) becomes 

r min =7 - 51xl0 ‘ 4 ^ CE p) 1/3 ( 27 ) 

The pusher diameter required to intercept the flux in a given cone angle 0 is ob- 
tained from the geometric relation (shown in sketch (c)). 



(28) 



Combining equations (26) to (28) yields the relation 


f - c 

d 2 1 

C 2 Vir 

d 2 + 2. 2 6X10” 6 (CE ) 2 / 2 

r 


(29) 


With the assumption that all the propellant mass has the mean velocity v, the upper 

limit to f c is 0. 5 in order to satisfy the momentum balance requirement. Equation (29), 

with the upper limit restriction of 0. 5 for f„, is shown plotted in figure 3 for four values 

o c m il l? 

of total pulse energy, E^ =4.18x10 , 4. 18xl0 A , 4. 18x10 1 , and 4. 18x10 joules (1, 

10, 100, and 1000 ton equivalents). (For reference, 1 gram of deuterium-tritium (DT) 
mixture fully reacted is equivalent to about 3.344x10** joules (80 tons TNT). ) The inter- 
relation among vehicle size (as reflected by pusher diameter), pulse energy size (Ep) and 
pulse unit design (as reflected by the collimation factor C) is evident. 

When the separation distance is maintained at the smallest value permitted by pres- 
sure limitations (the conditions imposed for fig. 3), one can note that there is a large 

improvement in the propellant interception factor f for (1) better collimation of the 

c 
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propellant, (2) smaller energy pulse value for the same degree of collimation, and (3) 
larger size vehicles (i. e. , larger pusher diameters), 

PROPELLANT-PUSHER INTERACTION 

When the flux of propellant arrives at the pusher plate, the initial, high velocity 
particles cause some sputtering. Some penetration into the material of the pusher also 
occurs. The effect of this initial bombardment is small compared to the eventual abla- 
tion caused by the arrival of the remainder of the pulse mass. The propellant that ar- 
rives is assumed to just lose its kinetic energy and form a hydrodynamic stagnation 
layer. The temperature of this stagnation layer is calculated by assuming that it reaches 
equilibrium through blackbody radiation back to the vacuum of space. 

Thus, when the energy flow rate relation (15) is used, the equilibrium can be ex- 
pressed as 


4 V AA 6 -(V *) 2 

oT* = — 2 — /_£) e r 

J/2.3 \ t / 

7 r ' r \ / 


(30) 


T s = 



11/4 


3/2 , 

.7 t ' err 


(31) 


A plot of this relation is shown in figure 4. Stagnation temperatures up to the 20- 
electron-volt range are typical. 

The formation of this high temperature layer causes the temperature of a pusher 
surface to rise quickly to the ablation level, The ablated gas then forms a protective 
layer and slows down subsequent ablation rates. 

At the temperature levels of the stagnation layer of gas, heat transfer to the pusher 
is mainly by radiation. Therefore, a surface material whose ablated products have a 
high absorptance for radiation in the frequency range characteristic of the temperature 
of the stagnation layer would be a distinct advantage in limiting ablation amounts. This 
situation is similar to that encountered in the gas-core nuclear rocket (refs. 10 and 11). 

In order to estimate the magnitude of the heat-transfer problem, the model of the 
interaction process shown in sketch (d) was assumed. The ablated vapor layer and the 
high-temperature layer from stagnating propellant were assumed to remain unmixed. 
The flow of heat through the ablated vapor layer was then looked at as a combined pro- 
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cess of radiation and conduction. It is assumed that the ablated vapor, of necessity, 
will be of material which has a high absorptance for radiant energy at the temperature 
of the stagnation layer. The one- dimensional equation describing the transfer of heat in 
an optically thick radiative- conductive medium with no heat source terms is 


_4o_8V +k 8?T 3T 

3a R 3x 2 C ax 2 p 3t 


(32) 


where a R is the adsorptance of vapor layer and k c the thermal conductivity of vapor 
layer. 

Numerical solutions to equation (32) may be obtained by assuming a boundary tem- 
perature history at x = 0 (given by the stagnation temperature relation (31)) and with the 
temperature gradient approaching zero as x becomes large. Typical temperature pro- 
files for these solutions are shown in figure 5. The rate of advance of the temperature 
"front" into the ablated vapor medium is also obtained. 

The inflection point of the temperature-distance relation occurs at about the temper- 
ature level for which equal rates of heat transfer occur by both radiation and conduction 
modes. This condition is given approximately by 


16oT 8 _ 
3a R 


(33) 


This relation becomes 


Tev ' 0 - 0128 ^/ 73 



(34) 


10 



when substituting for a and expressing temperature in electron volt units. The temper- 
ature level of the major portion of the ablated medium is below this inflection point tem- 
perature. This observation is noted here because of its relation to the properties of the 
ablated layer. 

The assumption that most of the radiant energy from the stagnation layer is absorbed 
in a relatively small thickness of the ablated vapor leads to the temperature profiles 
noted in figure 5, that is, a relatively steep temperature front that propagates through 
the medium. The similarity here to the gas-core rocket profiles is noted (refs. 10 and 
11). In the gas-core case, seeded gas is fed toward the advancing temperature front at 
such a rate as to maintain the steep temperature profile away from the wall of the rocket 
chamber. In the pusher case herein, ablated material is fed into the vapor state at a 
rate governed by the amount of heat that reaches the wall and the amount of energy re- 
quired to vaporize plate material and raise it to the temperature level of the ablated 
medium. 

In the following calculation of the amount of pusher ablation, it is assumed that heat 
transfer to the surface by radiation (though only a small fraction of the total radiant en- 
ergy available) is still the major mode of energy transfer. The incremental heat trans- 
fer per unit area, then, is 


dQ = aTg(t)e 


dt 


The increment of mass ablated per unit area is 


(35) 


dm Q = p dx 

a 



(36) 


Equations (35) and (36) were solved numerically for energy input using equation (30) 
when no ablated layer (x = 0) was present at t = 0 and when any decrease in x(t) result- 
ing from propagation of the temperature front into the ablated medium was neglected. 

The calculated ablation amounts will thus be lower than actual because of this latter as- 
sumption. 

Typical total ablation per unit area values are shown in figure 6. These were calcu- 

7 

lated using a value of 5x10 joules per kilogram for H„. The total ablation amount, how- 

a 

ever, is not especially sensitive to H„ as shown in figure 7. A factor of 10 change in 

<x 

H a causes only a 15 to 20 percent change in the total ablation amount under the conditions 
of these calculations. Figure 8 shows the range of energy absorption that goes into just 
ionization and thermal energy of the species for two metal vapors, iron (Fe) and uranium 
(U). Sublimation energy is, of course, also included in H . 

a 

The effect of separation distance on ablation rates is shown in figure 9 by a compari- 
son of curves. As the separation distance is increased, the energy arrival rate (per unit 
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area) decreases. However, the total interaction time increases. The net effect (fig. 9) 
of increasing the separation distance to four times the minimum separation distance 

12 

(eq. (28)) was to decrease the total ablation by only about 32 percent at E = 4. 18x10 

Q P 

joules (1000 tons) and 38 percent at Ep = 4.18xl0 joules (1 ton). 

The ablation calculations shown up to this point have been for conditions at the center 
of the pusher (/3 = 0)- Away from the center (/3 > 0), the ablation decreases as shown in 
figure 10. The same countering influence of the two factors, energy intensity and inter- 
action time, are present in the radial variation as in the separation distance variation 
(fig. 9). For Ep = 4. 18x10^ joules (1 ton), a 17 percent decrease in the ablation rate 
was calculated at 45° off the axis . 

Of special interest now is a calculation of the mass loss factor f m - Since the only 
such loss considered herein is that by means of ablation from the pusher surface, f m is 
given by equation (4). 

A cylindrical pusher plate of diameter d (sketch (e)) 



is now considered. The ratio of total ablation per pulse to propellant mass for the 
pusher plate of diameter d is 


5 

M„ 


r, fit, 

M„ / / H q 

p Jo Jo a 


0) -a R x 
e R dt 


dy 


P v «/Q 

The energy arrival rate at (z,/3) for a pulse with collimation factor C is 

3 


E t (z, /3) = Cp(z,/3) 


(37) 


( 38 ) 


Now, let t z = z/ v. Then t f = t z sec /3 and, using equation (16), 
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(39) 


E^z, 0 ) 


(_ % : 

- CE P V M e \t cos j3 y 


At this point the calculations are restricted to values of z corresponding to the maxi- 
mum pressure limitation (eq. (27)). When equation (27) is used, equation (39) becomes 


( y y 

E|(z, /3) = 4.23X10 8 v^ 6 e C ° S ^ 
Equation (37) is transformed to 


(40) 



by using the previous relations and putting the pulse energy in terms of Ep. 

The inner integral (with respect to t/t_) was solved numerically under the same as- 
sumptions as used in solving equations (35) and (36). The resulting values of the mass 
loss factor f for several combinations of variables are shown in figure 11. The 
curves stop at the limiting diameter where the propellant interception factor is 0. 5. 

The larger the pusher diameter, the smaller is the factor .f , other factors being 
constant. More pusher area is exposed to ablating conditions, while pulse mass Mp 
remains constant. The mass ablation loss factor is lower for smaller energy pulses for 
the same size pusher plate. This is because the separation distance z m is lower for 
lower energy pulses, and the pusher intercepts a greater fraction of the total energy flux 
from the pulse. 

The absorption coefficient a^/p has a marked effect on the mass loss factor since, 
as was noted before (fig. 6), the total ablation amount is nearly inversely proportional 
to a^/p for the conditions of these calculations. 


SPECIFIC IMPULSE 

A simplified model of the overall processes involved in the external nuclear pulse 
propulsion scheme has been assumed. Only the features of the propulsion system that 
affect the overall or effective specific impulse have been considered. The propellant 
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flow from the pulse energy source has been assumed to occur self-similarly from a point 
source. A collimation factor has been used to account for the concentration of more 
propellant into smaller flow cone angles than would occur simply from isotropic expan- 
sion. The hydrodynamic stagnation layer temperature is determined by the balance be- 
tween incoming kinetic energy and blackbody radiation back to space. Ablation of pusher 
surface material occurs by radiant heat transfer from the stagnation temperature through 
the ablated layer. A mean opacity for the ablated layer material is assumed. 

Using this model, two principal factors determining the effective specific impulse 
were calculated: the propellant loss factor f £ (eq. (29) and fig. 3), and the mass loss 
factor f m (eq. (4) and fig. 11). These two factors are combined to yield the ratio 

^sp^eff ~ *c ’ ^m ' ^sp^base 

Some typical plots of relation (42) are shown in figure 12 for a collimation factor of 3 and 
a R /p value of 10 square meters per kilogram. It is apparent in this figure that the op- 
timum specific impulse for smaller pusher diameters is attained with the small pulse 
energies. Under the conditions of figure 12 there is little or no improvement in effective 
specific impulse in going to higher mean propellant velocities at the smaller pulse energy 
levels. 


CONCLUDING REMARKS 

The calculation of effective specific impulse by the model used herein is felt to be 

optimistically high for a number of reasons. Most of the assumptions made relative to 

the calculation of the amount of ablation yield results believed to be lower than the true 

amount. Some examples are the following: 

7 

(1) The value of H„ used (5x10 J/kg) is at the high end of the expected range of 

a 

values (fig. 8). 

(2) The total heat transfer to the surface is greater than or equal to the radiant 
transfer portion assumed. 

(3) The decrease in px caused by the thermal wave propagation into the ablated 

-(a R /p)px 

medium is neglected in the exponential factor e 

(4) Mass loss other than the ablated amount M. may occur. 

The assumption of uniform distribution of propellant mass in the case of collimated 
flow leads to propellant interception factors that are higher than would be obtained for a 
more probable distribution. One would expect that in collimated flow mass would be 
more concentrated along the axis (/3 = 0). Thus, the stipulation of a maximum pressure 
limitation would mean that a greater separation distance would be required. This, then, 
would lead to lower f’ values for a given pulse energy and pusher diameter. However, 

V 
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the assumption that no impulse is contributed by the ablated vapor or by reexpansion of 
the stagnation layer leads to specific impulse values that are lower than the real case. 

The principal variables affecting the performance of the system are the following: 
Ep, the size of the energy pulse; r, the separation distance between the pulse and the 
pusher; v, the mean propellant velocity; C, the collimation factor; a^, the mean opacity 
of the ablated vapor; ( 3 , the angle from the pulse-pusher centerline; and H , the energy 

a 

absorbed during ablation. The effect of variations in these variables on the effective 
specific impulse of a given system is now discussed. The term "given system" is used 
here to mean one of constant diameter and a^/p value. 

(1) Pulse energy, Ep. There is an optimum pulse energy to yield the maximum ef- 
fective specific impulse with a given vehicle size (fig. 12). This result assumes that the 
separation distance is always kept at the minimum dictated by pressure limitation. The 
propellant interception factor increases (fig. 3) with decreased pulse energy. The total 
ablation amount decreases with smaller Ep (fig. 6); however, the ratio of the total abla- 
tion to the total propellant flow is greater at smaller Ep. Thus, the mass loss factor 

is lower at lower Ep values (fig. 11). 

(2) Separation distance, r. For a given pulse energy, locating the pusher farther 
from the pulse source cuts total ablation but not greatly (fig. 8). The decrease in energy 
arrival rate is counteracted by the increased interaction time. The big effect of r is, 
of course, a geometric one; that is, the greater the separation distance, the larger the 
required pusher diameter to intercept the same flux from the pulse. 

(3) Mean propellant velocity, v. There is no change in pressure with a change in v; 
the pressure is a function of the parameter CEp/r . However, the pulse energy re- 
quired for a given total impulse is directly proportional to v. 

There is no effect on total ablation with v in this model. While the energy arrival 
rate is proportional to v, the interaction time is inversely proportional to v; hence, 
the effects on total ablation exactly cancel. This result is obtained because the absorp- 
tance a^/p has been treated as a constant, independent of frequency. However, in a 
real system, the absorptance is a function of frequency (and hence the temperature of 
the radiation source). Some effect of v on ablation would be expected in such cases. 

(4) Collimation factor, C. The collimation factor, which has a marked effect on the 
propellant interception factor (fig. 3), is a problem associated with the pulse unit design. 
Thus, no estimate of the range of possible values is attempted here. A value of C = 3 
has been used herein for illustration purposes, but no implication about the probable at- 
tainment of this value is intended. 

(5) Rosseland mean opacity of ablated vapor, a^> This parameter is one of the most 
important ones affecting pusher ablation. The ratio of total ablated material to total 
propellant is nearly inversely proportional to a^/p under the conditions of interest 
herein. The importance of this property to the performance of the system makes the 
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selection of an ablating material and the correct determination of its a R of greatest 
impact. 

In the simplified approach of the model used, it was assumed that a R /p was con- 
stant. Experimental values of opacity of materials under the conditions of temperature 
and pressure encountered in the nuclear pulse interaction are not readily available for 
many materials. Some plots of the variation of the Rosseland mean opacity for several 
materials are shown in figures 13 and 14 (refs. 12 to 14). 

One can note by a comparison of figures 1 and 5 that most of the ablation has oc- 
curred considerably before the peak pressure (restricted to 6900 atm) is attained. The 
range of a R /p values indicated for materials in figures 13 and 14 would be generally 
below 10^ square meters per kilogram during most of the ablation process. 

The use of a mean opacity value independent of frequency (hence independent of the 
temperature of the radiating source) in calculating ablation rates undoubtedly leads to 
some inaccuracies in conclusions, as pointed out in item (3). The results of the analysis 
do, however, serve to establish the required range for the opacity at operating conditions 
in order that total ablation amounts will be kept within an acceptable range. Whether 
this range of a R values can be achieved with otherwise suitable pusher surface mate- 
rials is an area requiring much more detailed consideration. 

(6) Angular variation, 0. There is considerable variation in the flow properties 
(such as pressure, density, etc. ) from the center of a pusher to the edge (fig. 2). At 
the time the peak pressure occurs at the center (/3 = 0), the pressure at 30° off center is 
about 0. 55 times the peak value. The variation with /3 is greater than this at times 
prior to the occurrence of the peak value, and is less than this at times after the peak 
values have been reached. These time and radius variations in the pressure, for in- 
stance, undoubtedly complicate the pusher design to withstand the accompanying corre- 
sponding stresses encountered. 

(7) Energy absorbed during ablation, H a . An increase in H a decreases the total 

ablation as expected, but not nearly in direct proportion (fig. 6). At a pulse energy E„ 

9 P 

of 4. 18x10 joules (1 ton), for example, an increase in EL of a factor of 10 caused only 

cl 

a 15 to 20 percent decrease in total ablation. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, July 26, 1972, 

112-28. 
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APPENDIX - SYMBOLS 


Ross eland mean opacity, m” 
collimation factor 
heat capacity, J/(m)(sec)(K) 
pusher diameter, m 
energy, J 

pulse energy, J (1 ton TNT = 4. 18x10® J) 

2 

energy arrival rate, J/m sec 

2 

energy arrival rate at interface, x = 0 (see sketch (d)), J/(m )(sec) 

fraction of pulse mass that intercepts pusher 

ablation mass loss factor, defined in eq. (4) 

fraction of total mass flow that is included within cone angle, 6 

o 

gravitational acceleration, m/sec 
enthalpy, J/kg 

energy required to ablate pusher material, J/kg 

2 

impulse per unit area in z-direction at angle /3, N-sec/m 

total impulse in z-direction within included angle, 6 , N-sec 

specific impulse, sec 

base specific impulse, v/g, sec 

effective specific impulse, sec 

thermal conductivity, J/(m)(sec)(K) 

total mass ablated from pusher per pulse, kg 

propellant mass per pulse, kg 

o 

mass per unit area ablated from pusher, kg/m 

propellant flow rate per unit area, kg/(m )(sec) 

pressure, N/m^ (1 atm = 1.013xl0 5 N/m^) 

2 

maximum pressure, N/m 

2 

heat transferred per unit area, J/m 
radius, distance from pulse center, m 



r min minimum radius, m 

T temperature, K 

T g y stagnation temperature, eV 

T g stagnation temperature, K 

t time after pulse initiation, sec 

t m time at which maximum value occurs, sec 

t y reference time, r/ v, sec 

t z reference time, z/ v, sec 

Vj ionization potential, eV 

v velocity, m/sec 

v 2 z- component of velocity, m/sec 

v mean propellant velocity, m/sec 

W a weight of material ablated per pulse, N 

Wp weight of propellant per pulse, N 

w p propellant weight flow per unit area, kg/(m)(sec 2 3 ) 

x thickness of ablated layer, m 

y distance from 0 = 0 line in plane perpendicular to z, m 

z distance, pulse center to pusher plane, m 

z m minimum separation distance from pressure limitation, m 

0 angle from centerline 

6 included cone angle 

O 

p density, kg/m 

q 

P m maximum density, kg/m 

2 4 

a Stefan-Boltzmann constant, W/m K 
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Density at angle p 
Maximum density at p = 
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Figure 1. - Time variation of conditions at pusher location E p /r 3 of 2. 37xl0 9 joules per cubic 
meter. v 


Time, Time, 

t/t z t/i z 



(a) Density, p. (b) Pressure, p. 


Figure 2. - Radial variation of properties on pusher plate. 






Propellant interception factor, f ( 



Pusher diameter, m 


(c) Pulse energy, E p , 4. lSxlO 11 (d) Pulse energy, Ep, 4. 18xl0 12 
joules (100 tons). joules (1000 tons). 


Figure 3. - Propellant interception factor variation with pusher di- 
ameter, pulse energy, and collimation factor. 



0 2 4 

Time, t/t r 


Figure 4. - Time variation of stagnation temper- 
ature. 



0 .04 _.08 .12 

Ablated layer depth, px, kg/m 2 


Figure 5. - Temperature profiles in ablated medium. Pulse 
energy, Ep, 4. ISxlO 11 joules (100 tons); collimation fac- 
tor, C, 3; mean propellant velocity, v, 1. 5xl0 5 meters 
per second; absorption coefficient, a^/p, 10 3 square me- 
ters per kilogram; thermal conductivity, k c , 0.42 joule 
per meter per second per Kelvin degree. 


21 



Cumulative ablation, kg/m 


•02 1 


Pulse energy, 

Ep, 

J (tons) 



(a) Absorption coefficient, aR/p, 10^ square 
meters per kilogram. 




(b) Absorption coefficient, aR/p, 10 2 square (c) Absorption coefficient, aR/p, 10 square 

meters per kilogram. meters per kilogram. 

Figure 6. - Cumulative ablation per unit area as function of time for an energy required to ablate pusher 
material H a of 5xl0 7 joules per kilogram. Distance from pulse center equal to minimum radius, r = r m 
collimation factor, C, 3. 



(a) Absorption coefficient, aR/p, 10 3 square <b) Absorption coefficient, a R /p, 10 4 

meters per kilogram. square meters per kilog ram. 


Figure 7. - Cumulative ablation per unit area as function of time for pulse energy E p of 4. 18xl0 9 
joules (1 ton). Distance from pulse center equal to minimum radius, r - r m ; collimation factor, 
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Mass loss factor, 


Pulse energy, 

Ep. 

J (tons) 


Pulse energy, 
E P' 

J (tons) 


Pulse energy, 

Ep, 

J (tons) 




(a-1) Mean propellant velocity, 
v, IxlO 5 meters per second. 


(a-2) Mean propellant velocity, 
v, 1. 5x10* meters per second. 


(a-3) Mean propellant velocity, 
v, 2xl0 5 meters per second. 


(a) Absorption coefficient, a R /p, 10 2 square meters per kilogram. 



(b-1) Me^n propellant velocity, 
v, IxlO 5 meters per second. 


(b-2) Mean propellant velocity, 
v, 1. 5xl(P meters per second. 


(b-3) Mean propellant velocity, 
v, 2xl0 5 meters per second. 


(b) Absorption coefficient, a R /p, 10 3 square meters per kilogram. 


.4. 18xl0 10 (10) 


4. 18xI0 9 (1) 



(c-1) Mean propellant velocity, 
v, IxlO 5 meters per second. 



(c-2) Mean propellant velocity, 
v, 1. 5xl(P meters per second. 


I I I 

0 10 20 

(c-3) Mean propellant velocity, 
v, 2xl0 5 meters per second. 


(c) Absorption coefficient, 10 4 square meters per kilogram. 

Figure 11. - Mass loss factor variation with pusher diameter, pulse energy, mean propellant velocity and opacity. Collimation 
factor, C, 3. 
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Effective specific impulse, sec 


Pulse energy, 



0 10 20 0 10 20 

Pusher diameter, m 



(a) Base specific impulse, <b) Base specific impulse, (c) Base specific impulse, 

(Isp)base' 10 000 seconds. ^sp^ase* 15 000 seconds. (I spWe' 20 000 seconds. 

Figure 12. - Effective specific impulse variation with pusher diameter, pulse energy, and mean 
propellant velocity. Coilimation factor, C, 3; absorption coefficient, a R /p, 10 2 square meters 
per kilogram. 



Pressure, atm 


(a) Carbon. (b) Beryllium. (c) Hydrocarbon. 

Figure 13. -Variation of absorption parameter with pressure and temperature for 
several materia is. 
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Temperature, K 

Figure 14, - Variation of Rosseland mean 
opacity of uranium with temperature 
and pressure. 
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